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A new estimation of the interstellar proton spectrum 
is made in which the source term of primary protons 
is taken from shock acceleration theory and the 
cosmic-ray propagation calculation is based on our 
proposed nonuniform galactic disk model. 

It appers that above 10 GeV the interstellar proton spectrum can be 
determined with an absolute accuracy of ±20%, However, at lower energies 
the divergence among various spectral estimations is much more serious; 
this is due to the existence of solar modulation. Thus one needs to 
deduce the interstellar proton spectrum through his own demodulation 
calculation. The deduced proton spectrum is obviously dependent upon 
the model of solar modulation used. 

Recently,, an interesting attempt for deducing the interstellar 
proton spectrum is reported(l), in which a source term of primary protons 
suggested by shock acceleration theory is adopted and the cosmic-ray 
propagation calculation is based on the leaky box model. Since from our 
recent analysis of the high-energy electron spectrum(2) we are aware that 
the dominant part of observed protons comes from the dence H cloud region, 
it is argued that the corresponding propagation calculation should be 
based on the nonuniform galactic disk (NUGD) model(3). 

Here we will present a brief introduction to the NUGD model. 

According to it the observed cosmic rays contain two components: the 
distant component and the local component. As shown in Fig. 1 by the 
double-line arrows, the distant (left) component of cosmic rays, starting 
from the cosmic-ray confinement volume above the H^ cloud region (Box 2) 
can reach the solar neighbourhood by escaping into the H 2 cloud region 
(Box II), then by propagating along the magnetic tube (Box I) inside 
which the solar system is located. On the other hand, the local (right) 
component of cosmic rays originates from their confinement volume in the 
solar vicinity (Box 1). The fraction, € , of locally produced protons in 
the observed proton flux is found to be only 5±1 %(2). 

The main problem in extending the application range of the NUGD model 
to the low energy range is caused by a lack of knowledge of the convection 
velocity (V) of cosmic rays. Thus at low energies one is unable to 
determine a precise shape of the cosmic-ray escape pathlength (A ) in 
Box II or to estimate the cosmic-ray intensity variation along the 
magnetic tube (Box I) . Nevertheless, it is noted that if we limit 
ourselves to interstellar protons with kinetic energy, T ^ 1 GeV, to a 
great extent we can avoid the difficulties shown above. ^Actually, the 
proton intensity in the H 2 cloud region (Box II) should be the same as in 
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Box 1 or Box 2, because it is found experimentally that there is a lack 
of any cosmic-ray gradient, at least in the inner Galaxy (see Ref. 4). 
Note that at any given rigidity (R) protons are the dominant component of 
cosmic rays. Hence we do not need to use the unknown parameter to 

estimate the proton intensity in Box II (hereafter the subscripts 1, 

2, I and II represent the quantities referred to Boxes 1, 2, I and II 
respectively) . 

Regarding the A value in Box 1 or Box 2 we will use the empirical 
relationship( 1 ) deduced from measured data on heavy nuclei. Furthermore, 
in Box I the pathlength (x^.) of cosmic rays to reach the solar 
neighbourhood has been es timated( 2) . It is noticeable that even if 
cosmic-ray particles lack any convection motion in Box I, protons with 
T ^ 1 GeV should have a value of x le$s than 0.3 A , where A is the 
mian inelastic interaction length of interstellar prBtons, to ^ reach the 
solar neighbourhood. As a result, the proton intensity variation along 
the magnetic tube should be insignificant, and we will only consider two 
extreme cases (V= 0 and V= 300 km/s) to estimate the range of variation 
of the proton spectrum. 


need 


In order to deduce the proton intensity N ^ Box * or Box 2 we 
to solve the continuity equation of primary protons, 
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where q is the source term of primary protons and dN /dT = 1/T' is the 
energy ^distribution of protons after their inelastic^ ^ ^ 

interactions. From shock acceleration theory we have 


q pl2 = k 0p P 


-(2+7) 


( 2 ) 


where p is the proton momentum, k^ is a constant, = 0.05 and p c is the 
proton velocity. Since at low energies the power law approximat?on of 
N ^ an d the constant approximation of A 1 cannot be used, an iteration 
^ procedure is used to obtain the numerical solution of N ^ * 
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Fig. 2 

Further, the propagation of interstellar protons along the magnetic 
tube (Box I) can be described by using a slab model, 


(T ,x T ) N _(T ,x_) f oo 1 dN (T ,T' ) 
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where the initial value of N ^ should be N ^ after any possibly 
adiabatic deceleration in thE assumed bounSary layer s^ (see Fig. 1)0). 
At lower energies we also have 

* T - 4x . R _0 ' 7 ( l+(l+2q ) l, 2 r 2 , <4) 

x x u s 

where x IQ is a constant, q g = V^/Kj and is the dif f usion coefficient 
of cosmic rays. Similarly^ a numerical method is developed to obtain the 
solution of Eq. (3) in the solar neighbourhood, N g . Thus the proton 
intensity predicted for the NUGD model should be 

N = (1-e) N +£N , (5) 

PP ps pl2 

where t- 5±1 % as deduced from Ref. (2). 


It is found that the effect of taking various values of V is 
negligible, and the resultant cosmic-ray proton spectrum in the solar 
neighbourhood can be parametrized as, 7( . 
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C, = 15.0 - 6.05lnT + 2.841n T + 0.1691n T . 
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The predicted j spectrum is not only in accordance with the 
existing proton data P !t T ^ 10 GeV, but also is consistent with the low- 
energy proton data measured in the period of solar minimum(5) (see Fig. 2). 
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The modulated spectra in Fig. 2 are obtained from our predicted j 
spectrum by using the force field approximation. It appears that PP the 
measured data are consistent with our modulated spectrum with <t> - 200±200 
MV, where (p is the mean energy loss per nucleon in the heliosphere for 
cosmic-ray nuclei(6). 

Finally, in Fig. 3 our deduced interstellar proton spectrum (TA) is 
compared with the results of other authors (0P( 1 ) ,M0( 7) , ZU( 8) and JK(9)). 
It is noticeable that our proton speccrum is close to the demodulated 
spectra recently suggested in Refs. (8) and (9), indicating the general 
reasonableness of the demodulation processes used in these works. In 
conclusion, the parametrized expression of our deduced interstellar 
proton specrum should be useful in various cosmic-ray calculations, in 
particular in the estimation of the production rates of cosmic-ray 
secondary particles in interstellar spase. 
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